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A UNIFIED TWO-DIMENSIONAL APPROACH TO THE CALCULATION OF THREE-DIMENSIONAL 
HYPERSONIC FLOWS, WITH APPLICATION TO BODIES OF REVOLUTION ^ 

By A. J. EaoEEs, Je., and Raymond C. Savin 


SUMMARY 

A simplijied two-dimensional method jor calcidaMng three- 
dimensional steady and nonsteady hypersonic flows of an 
inoiscid (non-heat-condncting) gas is deduced from character- 
istics theory. This method is appropriately termed a generalized 
shock-expansion method. It is demonstrated that the method is 
applicable when disturbances associated with the divergence of 
streamlines in planes tangent to a surface are of secondary 
importance compared to those associated with the curvature of 
streamlines in planes normal to the surface. When this con- 
dition is met, surface streamlines may he treated as geodesics, 
which, in turn, may be related to the geometry of the surface. 

It is inquired further if the two-dimensionality of inviscid 
hypersonic flows has a counterpart in hypersonic boundary- 
layer flows. This question is answered in the affirmative, 
thereby permitting a unified two-dimensional approach to 
three-dimensional hypersonic flows. 

This concept is applied to bodies of revolution in steady Mght 
and, with the assumption that flow at the vertex is conical, 
approximate solutions for the flow field are obtained for values 
of the hypersonic similarity parameter (i. e., the ratio of the 
free-stream Mach number to the fluoness ratio of the body) 
greater than about 1 and for small angles of attack. Surface 
streamlines are approximated by meridian lines and the flow 
field is calculated in meridian planes. Simple explicit expres- 
sions are obtained for the surface Mach numbers and pressures 
in the special case of slender bodies. 

The validity of theory is checked by comparison with surface 
pressures and shock-wave shapes obtained experimentally at 
Mach numbers from S.OO to ■ 6. SO and angles of attack up to 
15° for two ogives having fineness radios of S and 5. At the lower 
angles of attack, theory and experiment approach agreement 
when the hypersonic similarity parameter is in the neighborhood 
of 1 or greater. At the larger angles of attack, theory tends to 
break down noticeably on the leeward sides of the bodies. 

INTRODUCTION 

The calculation of flows about objects, primarily missiles, 
traveling at high supersonic speeds is now generally accepted 
as a matter of more than academic interest. The difficulty 
of these calculations stems in lai^e part from the fact that 
at such high speeds disturbance velocities are not necessarily 
small compared to the velocity of sound, nor are entropy 
gradients necessarily negligible in the disturbed flow fleld 
about a body, even though it may be of normal slenderness. 
Thus, for example, the relatively simple linear theory, which 
has proven so valuable in studying flows at low supersonic 


speeds, loses much of its utihty in the study of bigh-super- 
sonic-speed flows. In the quest for methods especi^y suited 
to calculating high-supersonic-speed flows, notable progress 
has been made in the development of similarity laws relating 
the flows about slender three-dimensional shapes in both 
steady (see refs. 1, 2, and 3) and nonsteady motion (see 
refs. 4 and 5) . Steady two-dimensional flows have received 
perhaps the greatest attention from the standpoint of calcu- 
lating specific flow flelds, and it would seem that with tools 
ran^g from the characteristics method (see, e. g., refs. 6 
and 7) to the generalized shock-expansion method (ref. 7) 
the problem is reasonably well in hand, at least insofar as 
inviscid, continuum flow is concerned. A more or less 
analogous situation exists with regard to the nonlifting body 
of revolution (see, e. g., refs. 6, 8, 9, and 10) although it 
seems that only in the case of the cone has a method (ref. 10) 
of simplicity comparable to that of the linear theoiy been 
developed for calcxilating the whole flow field. 

When one departs from these relatively simple flows, the 
number of tools for carrying out practical calculations de- 
creases sharply. Thus, for example, in the category of 
inclined bodies of revolution, it appears that only bodies at 
small angles of attack have been handled adequately, usually 
by either the method of characteristics or some other step- 
by-step calculative procedure (see, e. g., refs. 6, and 11 
throu^ 14) . In the case of steady flow about genca^ three- 
dimensional shapes, aside from Newtonian flow concepts, 
which are strictly applicable at Mach numbers exceeding aU 
limits, only the characteristics method has apparently thus 
far received serious attention (refs. 15, 16, 17, and 18). It 
is true, of course, that the method is tedious and time con- 
suming to apply, but the relatively exact solutions obtained 
provide a valuable check against the predictions of more 
approximate but simpler theories. In addition, however, as 
demonstrated in reference 10, a study of the compatibility 
equations of the characteristics method can prove useful in 
determining simplified methods for calculating more complex 
flow fields. 

With these points in mind, it is first imdertaken in the 
present report to redevelop characteristics theory in a form 
which enables us to obtain a simplified two-dimensional 
method for calculating both steady and nonsteady hypersonic 
flows about three-dimensional shapes. Viscous flows are 
then considered and it is demonstrated that the two-dimen- 
sional character of inviscid hypersonic flows has a coimterpart 
in hjrpersonic boundary-layer flows. The validity of the 
analytical methods of this paper is checked by comparing 


1 Supersedes NAOA TN 2811 entitled "On tire Calcnlation of Jlcm- Abont Objects TnsTellng at High Supersonic Speeds," by A. J. Eggera, It., 18S2. 
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the predictions of theory with experimental results for the 
surface pressures and how shock waves of lifting and non- 
lifting bodies of revolution at Mach numbers from 3.00 to 6.30. 
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NOTATION 


local speed of sound 
axial-force coefficient, 


axial force 


•<l) 


normal-force coefficient, 


normal force 


2-^ I j 


pitching-moment coefficient, 

moment about body vertex 



specific heat at constant pressure 
specific heat at constant volume 
characteristic coordinates in X-Z plane (<7i, is 
positively inclined with respect to X) 

pressure coefficient, - — — 


maximum diameter of body of revolution 
hypersonic similarity parameter, j 


characteristic body length (measmed from vertex 
to most forward point of maximum diameter) 
Mach number (ratio of local velocity to local speed 
of sound) 
static pressure 
total presstue 

free-stream dynamic pressure 

entropy 

time 

components of fluid velocity along the X, Y, and Z 
axes, respectively 

rectangular coordinates along the X, Y, and Z axes, 
respectively 
cylindrical coordinates 

center-of-pressure position (measmred from body 
vertex) 

angle of attack 

ratio of specific heats, ^ 

angle between X axis and tangent to projection of 
streamline (or pathline) in X-Z plane 
semivertex angle of body 

angle between X axis and tangent to projection of 
streamline (or pathline) in X-Y plane 

Mach angle, sin“* ^ 

ray angle for Prandtl-Meyer flow 
mass density 


SUBSCRIPTS 

free-stream conditions 

jconditions at different points in the flow fleld 
conditions on the surface at the vertex of a body 


s conditions immediately behind the shock wave at 

the vertex of a body 

mVISCID FLOVT 


This study proceeds from the Euler momentum equations, 

( 1 ) 


dii , bu . bv, , bu 1 bp 


dp , dp , dp , dp 1 dp 
bt bx~^^ by~^^ dz^ p by 

bw , bw , bw , bw 1 bp 


the continuity equation. 


dp I d (pa) d (pp) ■ d (pw) 
bt'^ bx by bz ^ 


the equation of state. 


P=p(p,S) 


and the energy equation. 


bS , bS , bS , dS_„ 


(2) 

(3) 

(4) 
( 6 ) 
( 6 ) 


where u, v, and id are the components of velocity at time t 
along the X, Y, and Z axes, respectively, of an element of the 
fluid of density p, static pressme p, and entropy S.* To put 
these expressions in a more tractable form, it is convenient 
to aline the X axis at time t with the direction of the resultant 
velocity at the origin of the coordinate system. Thus 
equations (1) through (4) and equation (6) simplify, respec- 
tively, in the region of the origin, to 


bu , dw , 1 dp „ 
d^ bx pbx 

(7) 

bv , dp , 1 dp . 

■3^+'*^ 3T+- 
bt bx p by 

(8) 

bw , bw ,1 bp „ 
-^^Hru — h- -^=0 
bt bx p bz 

(9) 

bp , bp , /bu , bv , bw\ „ 

(10) 

as, as „ 

ai) 


which relations are basic to the subsequent analysis. 

- STEADY FLOW 

Characteristics theory. — Compatibility relations describing 
the behavior of fluid properties along characteristic lines 

* For certain calcnlatioiu it may bo desirable to proceed Irom moro gonerol equations which 
indodo etldrts of heat and mass addition to (or subtraction from) the flow os well os olTcots of 
impressed forces (e.g.> gravitational or magnetlo). Such a procedure may easily bo developed 
from that presented here by following the method of Qudorloy (roferonco 0) for two> 
dfmpTHrifmfll flow. 
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in supersonic flow may, f course, be obtained by proceeding 
formally with the theory of characteristics for the quasi- 
linear partial differential equations which depict the flow. 
In the interests of simplifying both the derivation of these 
relations and their resultant forms, however, it seems desir- 
able to proceed in a more intuitive manner, assuming a 
priori that the pertinent characteristic lines are Mach lines, 
and utflizing the implication from two-dimensional flow 
studies that perhaps the most convenient dependent vari- 
ables are pressure and flow inclination angles. 

Now it is clear that in the case of steady flow all derivatives 
with respect to time disappear from the above relations. 
Thus, assuming there are no shock waves present in the 
region of the origin,® we may write, with the aid of equations 
(6) and (11), 

bp 5p 

bx~i>p 


op op^ 

5a: a* 5a: 


( 12 ) 


where a is the local speed of sound in the fluid. Combining 
equations (7), (10), and (12), there is then obtained the 
relation 


5p_ —pit® ri /5o 5 m)\~| 


(13) 


or, defining A as the angle between the X axis and the 
tangent to the projection of a streamline in the X-T plane, 
and, in an analogous manner, the angle 5 in the X-Z plane 
(see fig. 1), we have 


Streomline 



PiauBB 1. — Streamline projections in X-Y and X-Z planes. 


5p —pu^ /5A ■ 5g\ 
1 \5y”^5s/ 


(14) 


Transforming the derivatives with respect to x and z to 
derivatives in the characteristic or Ci, and C 2 , directions 
in the X-Z plane (Ci, is positively inclined with respect 
to X, thus bQ/^x=[M/{2^/M^)][^0/^Ol.+^0/^>C2,] and 
bQldz=(Mf2)[bQlbCu—^0I^Ou]) there results from this 
equation 


5p &p _ , pit® r 55 n 


In an analogous manner, there is obtained from equation (9) 
the relation 


dp dp — ptt® / 55 55 \ 

5(7u“56'2,~y#=:i 


(16) 


Adding these two expressions then yields 

5p _ — pu* r 55 . 1 /2>A\"] 

wMle subtracting yields 

dp _ pa® r 55 ^ /5A\~| 

1 AfW/J 


(17) 


(18) 


Equations (17) and (18) are compatibility equations for 
characteristic or Mach lines in the X-Z plane.'* Indeed, if 
it is further required that the X-Z plane be the osculating 
plane of the streamline passing through the origin, that is, 
the plane containing the principal radius of curvature and 
tangent to this streamline (at the origin), then these equa- 
tions are the essential relations for determining pressure 
and flow inclination throughout a flow field. This point 
becomes evident when it is observed that, with the imposed 
requirement (viz, 5A/5r=0), the additional information 
derived from studying flow in the X-Y plane is simply 
that deduced from equation (8), or, as would be expected, 



(19) 


In order to construct a flow field, however, it is necessary 
to know the manner in which the osculating plane rotates 
and, correspondingly, how the principal curvature varies as 
we proceed along a streamline. This information is obtained 
from equations (2) and (3). Differentiation with respect 
to X yields 


5®A ^ ^ /5p\ 5A 55 

5^~ pu® 5a: \dp/ dz dx 


( 20 ) 


and 


5®5_ 1 5 / 55 . 3 5u. ■ 1 5p\ 55 

5a:®~ pTi® 5a: \5z/ \dz^H dx po!‘ dx) 5a: ^ 


respectively. 

These and the previously derived expressions foim the 
basis of a characteristics theory for steady three-dimensional 
flows (see ref. 20). Consistent with the objectives of this 
paper, however, we are interested in these results as they lead 
us to a more approximate but, by the same token, a more 
simplified method of calculating the three-dimensional flow 
of a gas at high supersonic speeds. 

Simplified two-dimensional theory. — ^It is well at the out- 
set of this analysis to establish, insofar as is practicable, the 
type of flows to be treated. In this connection, it is con- 
venient to employ the hypersonic similarity parameter 

e., the product of the flight Mach number and the thickness 
ratio of a body) as a measuring stick. In flows charac- 
terized by values of the hypersonic similarity parameter 
small compared to 1, that is, flows in which the body is 
extremely slender and hes close to the axis of the Mach cone, 
there is no apparent reason to believe that the linear theoiy 
will not be as useful an approximate method of calculation. 


• IX shook wnvcs are prejont, the apprepriate oblique shock equations are employed. 


< It is noted that these expressions contain not only derlTatives In the characteristio direc- 
tions hut also deriratlves with respect to the Independent variahle y. This type oi result is 
to be expected as pointed out by Oobum Crat 19). 
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as at low supersonic speeds. In flows characterized by values 
of the parameter up to about 1, the second-order theory first 
enunciated by Busemann (ref. 21) for airfoils and more 
recently generalized to three-dimensional flows by Van Dyke 
(ref. 9) and Moore (ref. 22) shoidd prove a useful approxima- 
tion. On the other hand, for flows aBout’ more or less 
arbitrary shapes, there is apparently no approximate method 
of calculation generally applicable with engineering accuracy 
at values of the hypersonic similarity parameter appreciably 
greater than 1. 

. In the limiting case of indefinitely high free-stream Mach 
number (and hence similarity parameter) and a ratio of 
specific heats equal to 1, we have the Newtonian impact 
theory (ref. 23) and its refined counteiqjart, accounting for 
centrifugal forces in the disturbed flow, developed fihret by 
Busemann (ref. 24) and more recently treated by Ivey, 
Elunker, and Bowen (ref. 25). The impact theory has been 
employed with some success by Grimminger, Williams, and 
Young (ref. 26) and others to predict smface pressmes on 
bodies of revolution at values of the similarity parameter 
appreciably greater than 1, although it should be remarked 
in passing that this success is in part, at least, fortuitous, as 
perhaps is best evidenced by the fact that the more exact 
theoiy (within the framework of the imderlying assumption 
of , y -^1) of Busemann is considerably less accurate 

under corresponding circmnstances. As shown in reference 7, 
neither the Newtonian impact nor the Busemann theory 
apply with good accuracy to airfoils except at values of the 
similarity parameter quite large compared to 1, correspond- 
ing, for example, in the case of thin airfoils to flight ^eeds 
considerably in excess of the escape speed at sea level. 
Perhaps the foremost shortcoming of these theories is, how- 
ever, that, irrespective of the shape to' which they are applied, 
they provide no information on the structure ® of the dis- 
turbed flow field which is, of comse, of finite extent adjacent 
to the surface at^flight Mach numbers presently of interest 
(say Mach numbers less than the escape Mach number at 
sea level). Such information is, for example, important to 
the determination of the flow about control surfaces and the 
like which may be located in this field. 

In view of the preceding discussion, it seems dear that in 
the high-supersonic-speed flight regime, a need for an 
approximate method of analysis hes in the realm of flows 
characterized by values of the hypersonic similarity param- 
eter greater than 1. An attempt will therefore be made to 
obtain a method meeting part of this need, attention being 
focused primarily on flows characterized by large values of 
the similarity parameter. To this end, it is convenient first 
to employ equation (14) rewritten in the form 


where 


pa* pi— 1 /dA~\~| 

1 V1+-D J -^W—l 

-r, bS/bCit 


( 22 ) 

(23) 


Now consider for the moment a surface streamline aimed 
in the x direction, and impose the requirement that the X — Z 

■ This consequeilce Is traceable primarfly to the asomption of 7*1 which leads to the 
well-known result that the distorbed flow field Is confined to an region adjacent 

to the stuface of a body.' 


plane be tangent to this streamline and normal to the sur- 
face at the point of tangency (the origin). The X — Y plane 
is then, of course, tangent to the surface at this point. Ob- 
serving the last term in the brackets on the right-hand side 
of equation (22), it is noted (see fig. 2) that 


Y 



T 



FiauBB 2. — Divergence of streamlines in tangent piano. 


dA 1 
dy r' 


where ri is the radius of curvature of the line normal to the 
projections of streamlines in the X-Y plane, and passing 
through the origin. At the high Mach nmnbers imder 
consideration, the disturbed flow field is confined to a region 
of small extent normal to the surface of a body; hence it 
may be expected that r' wUl be primarily a function of body 
shape and attitude. * Thi s being the case, it foUoAvs then 
that the term (l/vAP— l)(l/r'') will decrease in absolute 
magnitude with increasing Mach number of the flow about 
the body. Consider now the term (5S/&z)(l— D,/(l-fP,). 
We note that bSfdx=llR where B is the radius of curvature 
of the projection of a streamline in the X~Z plane and, 
by reasoning analogous to that used in considering r', is not 
expected to vary significantly with Mach number in the 
disturbed flow field. Let us assume for the moment that 
the quantity (1— 2?,)/(l-l-D,) is also relatively independent 
of Mach munber. With this assumption, it is clear that 
equation (22) approaches the equation for two-dimensional 
flow as the free-stream Mach number, and hence the hyper- 
sonic similarity parameter of the flow becomes large compared 
to 1. The compatibihty equations (eqs. (17) and (18)) ari 
affected in a similar manner; thus it is apparent that the 
flow when viewed in the X—Z plane approaches the two- 
dimensional type. In this case, however, as shown in refer 
ence 7, so long as the Mach munber and ratio of specific 
heats of the disturbed fluid are not too close to 1, D% is smal 
compared to 1, and hence the flow approaches the generalizec 
Prandtl-Meyer type (i. e., flow in which pressiuo anc 
inclination angle are approximately constant along emwee 
first-family Mach lines). Our flow equation may then b( 
written 


pnj? /dS\ 


(24 


where it is required exphcitly that 


Idx 


» 


1 

&A 

-^/AP-l 

dy 


(26 


° It b Inteiestlng to note tbat in ideal gas flows, r' becomes Justolnnetlon of these vorlabli 
as the valne of the hypeisonlo shniiarlty patameter becomes large compared to 1 (see wor 
of Oswatitseb, ref. 27, noting that hb results can readfly be extended to thrcOHlImonalom 
ideal gas flows nsing the obaraoteristics method of tbb paper). 
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or, in effect, that disturbances associated with the divergence 
of streamlines in tangent planes must he of secondary im- 
portance compared to those associated with the curvature of 
streamlines in planes normal to the surface. Upon closer 
o.vamination it can easily be deduced that this requirement 
stems directly from the continuity equation and the condition 
that it be of the two-dimensional type in the X-Z plane (i. e., 

in eq. (10)). 



» 

dv 

dz 


From t 


these considerations it appears that the conclusion 
of reference 10 that inviscid flow along streamlines down- 
stream of the nose of noninclined bodies of revolution travel- 
ing at high supersonic speeds may be of the Prandtl-Meyer 
type (in regions free of shodk waves) applies also to other 
steady three-dimensional flows. It is true, too, that in the 
latter case, just as in the former case, this conclusion is 
consistent with the predictions of the hypersonic similarity 
law for steady flow about slender shapes. 

One question remains to he considered, namely, where do 
the streamlines go in the disturbed flow? To clarify this 
matter, it is convenient to study further the implications of 
equation (25). For this purpose we combine equation (25) 
with the transformation equation 

dA M / Z) A dA 


dy 


__M / dA 

~~ 2 \dOif dOig) 


to obtain the relation 


>- 


M 

&A 

5A 

2VM®-1 

dOiy 

dOts 


(26) 


From this relation we deduce either that to the order of a 

da: 


number (curvature) small compared to ^ ^ 


or that 


and 


dA dA 

dOiy d(7jy 


dS 

.. M 

5A 


dx 

"2VM®-1 


- 


» ^ 

dA 


5a: 

2VM**-! 

dOi, 



(27) 


(28) 


Equation (27) implies vortical flow, however, which type of 
flow cannot he treated by the present analysis since equation 
(25) is violated.^ Equation (28) is then the requirement con- 
sistent with the basic assumptions of this analysis. Com- 
paring the relations of equation (28) with the transformation 
equation 


M- 

- ^ 1 



5a; 

2VM^-i' 


^dOi, 


leads one to the conclusion, however, that 


da;] 




dA 

da: 


(29) 


1 This conclusion Is particularly evident In the case of pure vortic^ flow, cff say vortical . 
flow with a superimposed uniform stream direct^ along the axis of the vortex, In which 
bp .... 

cases ^ eo, and hence equation (24) certainly does not follow from equation (231). 


or, in effect, that consistent with equation (25) geodesic lines 
can be treated as surface streamlines. With this informa- 
tion we are enabled to construct the flow field about a body, 
having once determined, for example, the flow in the region 
of the leading edge (or edges) thereof. This result follows 
since a geodesic line, and hence a streamline, on the surface 
is fixed, provided its direction at any point is given (see, e.g., 
ref. 28).® With this knowledge of the location of surface 
streamlines, flow in the planes tangent thereto and normal 
to the surface may be calculated approximately in the rela- 
tively thin region between the surface and hoimding shock 
waves, using the generalized shock-expansion method in the 
manner described in reference 7. 

A partial check on these observations is afforded by study- 
ing the flow about a swept airfoil. In this case flow at the 
surface may be calculated with good accuracy, using the 
shock-expansion method in combination with simple-sweep 
theory. For thin airfoils (on the surfaces of which the 
appropriate geodesics have essentially the direction of the 
free stream) the generalized shock-expansion method of this 
paper reduces to the slender-airfoil method of reference 7. 
Thus, in this case, it is evident from the results of reference 7 
that the generalized method will predict surface pressure 
coefl&cients in error by less than 10 percent, providing the 
component of free-stream Mach number normal to the 
leading edge is greater than about 3. It is of interest also 
to consider a ttdek airfoil to ascertain the accuracy with 
which this method applies to flow with appreciable curvature. 
To this end, surface pressure coefficients and streamlines 
have been calcrflated for a 20-percent-thick biconvex 
airfoil (at zero incidence) swept 60° and operating at Mach 
numbers of 10 and infinity ( 7 = 1 . 4 ). Conditions at the 
leading edge were determined from exact shock-wave rela- 
tions for both methods. The results of these calculations 
are presented in figure 3, and it is observed that the pressure 
distributions determined with the shock-expansion method 
for swept airfoils and the generalized shock-expansion method 
are in reasonably good agreement at both Mach numbers. 
The streamlines are also in reasonably good agreement over 
the forward portion of the airfoil, although, as would he 
expected, somewhat poorer results are obtained over the 
afterportion. It is not surprising, in view of the underlying 
assumptions of the generalized shock-expansion method, 
that it is generally more accurate at the highest Mach 
niunber. 

In the preceding discussion circumstances were deduced 
imder which steady flow at high supersonic speeds about 
three-dimensional shapes could he constructed approxi- 
mately, using the basic tools of two-dimensional supersonic 
flow analysis, namefy, the oblique shock equations of 
Rankine and Hugoniot and the comer expansion equations 
of Prandtl and Meyer. Several possible exceptions to these 
circumstances immediately come to mind. These include 
conical-type flows and flow in the region of the tip of a wing, 
or at the discontinuous jimcture of a wing and body, to 
mention a few. In such flows equation (25) may not he 
satisfied, in which case two-dimensional flow in planes 

• Ha sudden change of surface slope causes an oblliine shock wave or a concentrated Prandtl- 
Meyer type eipansJoli fan, the streamlines In the downstream direction are defined on tha 
basis of their flow direotlan Immediatoly following the discontinnity in slope. 
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normal to a surface cannot be expected.® It might be 
reasoned, therefore, that these flows cannot, in general, be 
treated by the proposed method. This observation may be 
correct; however, in one case investigated thus far in this 
connection, namely, flow in the region of the nose of non- 
lifting bodies of revolution (see ref. 10), it was found that 
although equation (25) is not satisfled, flow along stre amlin es 
is nevertheless of approximately the Prandtl-Meyer type. 
Thus we are led to expect that perhaps a less restrictive 
requirement than the satisfying of equation (26) may be 
imposed to insure that flow along streamlines is of this type. 
Such a requirement is in fact easily obtained by reconsidering 
equation (22) in the form 


< One may note that In some cases of this nature, tlie flow in oscolating planes of the stream* 
may be of the two-dimensional or even the aimpler Prandtl-Meyer type, although these 
planes may not be normal to the sorface. 



pU* 

“55 M 

/ 55 


5a; 

VM^-1 

L5a; VMS— 1 

Km. 

M by)] 


thus yielding 



M 55 


bx 

Vm® — 1 bO\t 

'M by 


(30) 

(31) 


It is evident that equation (31) embraces equation (26) as a 
special case and that Prandtl-Meyer flow obtams along 
streamlines if 


b5 1 &A 
Mby 


(32) 


to the order of a number small compared to 




~JT 


&a: 


This result imphes that although flow inclination angles are 
not necessarily constant along Ou lines, pressure is approxi- 
mately constant (see eq. (17)). 



Sponwise station, y/c 


Chordwise station, x/c 


(a) M„ = 10 

(b) iW„= => 

FianEB 3. — Comparison, of surface streamlines and pressure distributions calculated with the generalized shock-expansion method and 
the shock-expansion method for swept airfoils (biconvex airfoil section, thickness ratio=0.2, sweep angle=»60“). 
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It is clear that the increased generality of the above result 
has been obtained at some expense in our knowledge of the 
streamline flow pattern. For example, it is not now indi- 
cated that (within the framework of this analysis) surface 
streamlines may generally be taken as geodesics — additional 
knowledge of the flow must be had in order to determine 
these streamlines. If they are known, however, the calcula- 
tion of the whole flow field is materially facflitated by the 
above considerations. 

Thus far only steady flows have been considered. The 
next problem is to extend these considerations to nonsteady 
flows and some aspects of this matter will now be discussed. 

NONSTEADY FLOW 

The methods of analysis in this case are entirely analogous 
to those employed in the study of steady flow, the singular 
contrasting feature being that derivatives with respect to 
time in equations (1) through (11) cannot now be neglected. 
With this point in mind, only pertinent results are discussed 
below. 

Characteristics theory. — The compatibility equations relat- 
ing fluid properties along Mach lines may be written as 
follows: 


bp _ ~pu^ ( 
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The definition of the X-Z plane as the osculating plane of a 
pathline (streamline in steady flow) remains as before, hence 
equation (19) still applies in the X-Y plane in the region of 
the origin. The rotation of the osculating plane and varia- 
tion of the principal curvature of a pathline with motion 
along it are now, however, obtained with the aid of the 
relations 
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These equations are basic to characteristics theory in its 
application to three-dimensional nonsteady flows (see ref. 
20). However, as m the case of steady flows, they can best 
be employed for our purposes to obtain a simplified method 
of calculation. 


Simplified two-dimensional theory. — ^It is recalled that the 
essential simplification in our treatment of three-dimensional 
steady flows derived from the fact that they often appear 
locally two-dimensional and hence they can be treated with 
the generalized shock-expansion method of reference 7. In 
the foUowing study of nonsteady flows we will profit from this 
experience by anticipating that the desired simplified theory 
is again this shock-expansion method. Accordingly, our 
problem is reduced to that of detennining the conditions 
imder which the method can be applied to the calculation of 
nonsteady hypersonic flows. 

One condition is, for all practical purposes, self-evident; 
namely, the local Mach number of the disturbed flow must 
be everywhere large compared to 1. This requirement must 
manifest itself since, otherwise, nonsteady disturbances 
created an appreciable distance upstream and/or downstream 
of a particle could significantly influence its behavior in the 
disturbed flow field (see fig. 4, noting that in case of thick 
body, particle b is influenced by disturbances originating in 
particles a and c) and this situation would preclude the pos- 
sibility of Pi-andtl-Meyer type flow along pathlines. It fol- 
lows then that the shock-expansion method be applied 
only to thin or slender shapes (i. e., shapes producing flow 
deflections small compared to 1) at hypersonic speeds. 


Fluid pot tides — 



- Wove fronts of 

disturbonces \\ 

generated in portides \ X 



Shock wove 


M »l- 


Figiteb 4. — Propagation of nonsteady disturbances in flow about thick 
• and thin bodies. 


With this requirement in mind it is convenient to rewrite 
the compatibflity equations (33) and (34) in the form they 
assume when M'^1. Thus with slight rearranging we have 

bp ■ 1 —pu^ / 55 ■ 1 dS ■ 1 &A 1 bv\ 
bOu~^^ bt ^ M VStC'ii « bt^M by Mv? bt ) 

and 

dp . 1 £)p pu,^ / £>S ■ 1 dg ^ dA . 1 bu\ 

bOu 'u, bt \dC7i, ubt M by'Mu^ bt ) 

Now, (insistent with the requirement M^l, the term 
(1/iliu*) (da/dt) on the right in these equations can be 
neglected by comparison to the other terms. If, in addition, 
we define (after eq. (37)) the derivatives 

d d 1 

dOj^bOu u \dty 

and 

d d 1 /d\ 
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and note that now 
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then these equations can be combined to yield the pressure 
gradient along a pathline in the following form 




r ds/do^-] 
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(38) 


But is analagous to D, in steady flow; that is, in the 

event the term (1/il/) (jbA/dy) is negligible in this expression, 

identified with disturbances reflected from 

shock waves in the flow. Just as in the case of steady flow 
(see ref. 7), however, these reflected disturbances are of very 
small strength by comparison to the incident disturbances 

when M^l and so must he small compared to 1. 


Provided, then, that 


dSjdOii 


M»1 



dS 

„ 1 

bA 

> 

dx 


by 



equation (38) may be written “ 


(39) 
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dx~ M dx 
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which result implies, of course, Prandtl-Meyer flow along 
pathlines. It follows that equations (39) are sufficient con- 
ditions under which the generalized shock-expansion method 
can be used to c-alculate nonsteady hypersonic flows. 

When these conditions are satisfled, we note, by analogy 
to the steady flow case, that pathlines in the surfaces swept 
out by elements of fluid adjacent to shapes in 'nonsteady 
motion are approximated by geodesics or, even simpler, lines 
of cmwature of these surfaces. It is not to he imphed, of 
course, that pathlines must always be such curves in order 
for fluid properties to behave as in Prandtl-Meyer flow. In 
fact, again, just as in the case of steady flow, if the condition 


dS 

dx 




dS 

■ rdA 

dCu 

by 


is satisfied rather than the second 


of 


equations (39), pathlines are not necessarily geodesics (or 
lines of curvature) eveii' though the first' of equations (39) 
and hence equation (40) holds along these lines. 

One notes that within the framework of this approximate 
analysis, the calculation of nonsteady flows at the surface 
of slender bodies traveling at high supersonic speeds should 
no t pro ve unduly difficult. To illustrate, consider an oscfllating 
airfoil as shown in figure 5. The presstire at any point 
along the pathline shown is readily deduced by simply 


integrating equation (40) along this line from the leading 
edge of the airfoil to the point in question. The whole 
flow field as a fxmction of time may he calculated by 
employing the generalized shock-expansion method for 
steady flows (see ref. 7) in a series of planes located small 
distances apart in time. This example serves to emphasize 
that, in general, the time history of fluid elements must he 
known, at least to the extent of fixing their initial flow direc- 
tion and entropy.” It is also evident that again, as in the 
case of steady flow, the general results of the analysis are 
consistent with the predictions of the hypersonic similarity 
law for nonsteady flows about slender related shapes (ref. 6). 

These considerations complete our general treatment of 
inviscid hypersonic flows. It is appropriate to turn next to 
eflfects of viscosity as they relate to the hypersonic boundary 
layer. 

VISCOUS FLOW— THE HYPERSONIC BOUNDARY LAYER 

The arguments presented here are concerned Avith the 
steady hypersonic boundary layer, and they wfll he, for the 
most part, physical.^ Furthermore, they wall appear as 
natural extensions of conclusions reached in our study of 
inviscid hypersonic flow. Let us reconsider, then, the motion 
of the inviscid fluid. We have established that this motion 
is, under certain well-defined circumstances, confined locally 
to planes normal to the surface of a body and tangent to 
surface streamlines. Correspondingly,, there is no sensible 
momentum transfer across these planes. Now if viscous 
forces are set up in the flow hounding the surface, we recog- 
nize that they wfll act to resist the motion of the fluid — that 
is, the motion in the normal planM. Evidently, then, these 
forces act in the same planes of local two-dimensional flow 
as the pressure forces, and it must follow, of course, that 
resultant changes in momentxim of the fluid also occur in 
these planes. 

Consider now the changes in energy of the disturbed fluid. 
These changes can he brought about by viscous or dissipa- 
tive work, pressure work, heat convection, and heat conduc- 
tion.^® It was just found, however, that the forces doing 
work act in the normal planes; hence we conclude that the 
corresponding changes in energy occm in these planes. 
Similarly, heat is convected in the normal planes, since mass 
is convected in these planes. Finally, we conclude also that 
heat is conducted locally; in the normal planes inasmuch ns 
the temperature gradients set up by the action of viscous 
forces are confined primarily to these planes.” Evidently, 
then, changes in energy of the fluid can be treated locally as a 
two-dimensional phenomenon in planes normal to the smface 
of a body. 


Airfoil of time /g 


PothTme of porficle striking 
leading edge at time /. 


Airfoil at time / 



- — Shock wave 


.Direction of rotation 
of airfoil 


Figuee 6. — Oscillating airfoil in hypersonic flow. 

® The derivation of equation (40j presented here has tlie v irtu e of being both more complete 
and more general than the corresponding derivation of reference 20. 


11 In the special case of slender airfoils for which the hypersonic similarity parameter of the 
Sow is less than 1, entropy gradients in the disturbed Sow can frequently be negleetod, with 
the advantage of relaxing these two conditions and thus substantially simplifying the prohlom 
(see, e. g., ret 29). 

“ Although not presented, corresponding mathematical arguments have hcon pursued 
using the Navler.Stokes and energy equations, and the Onal results confirm those obtained 
here. It is indicated, too, that these results may apply also to nonsteady boundary-layer 
Bows. 

i» Radiation and absorption may, of course, also contribute to the energy changes: however, 
It is beyond the scope of this paper to consider these phenomena. 

X One might conceive of severe temperature gradients being imposed at tho wall boundary 
by, for example, extremely nonunlform surface cooling. Such gradients, If transverse to 
streamlines, would naturally Invalidate this argument 
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Thus far we have been concerned mainly with forces and 
their relation to the momentum and energy of the fluid. 
The question of conserving mass remains to be investigated. 
It will be recalled that the requirement of conservation of 
mass was the essential factor which determined when the 
generalized shock-expansion method could be employed to 
calculate three-dimensional flows. This requirement is 
physically (and mathematically) the same, independent of 
whether or not viscous forces come into play. We conclude 
then that for the purposes of this study, equation (25) can be 
used to determine when the three-dimensional boundary layer 
can be calculated with two-dimensional equations. From 
equation (25) it is indicated that the boimdary layer must be 
largely hypersonic if this calculation is to be permissible. 
It is not to be imphed, however, that the boxmdary layer 
always becomes two-dimensional, as on an airfoil, if the 
stream Mach number is mode extremely large. For example, 
in the case of axial flow about tbe right circular cone, equa- 
tion (25) is violated independent of Mach number (just as 
with inviscid flow) and we must use something like the 
Mangier transformation (ref. 30) in the boundary-layer 
calculations. On the other hand, if the body, instead of 
being conical, is curved in the stream direction, then it is 
indicated that the boundary-layer flow should approach the 
two-dimensional type with increasing Mach number. 

This discussion completes our arguments regarding the 
two-dimensionality of three-dimensional hypersonic flows. 
Attention is turned next to a practical application of this 
concept. 

APPLICATION OF THEORY TO BODIES OP REVOLUTION IN 
STEADY FLIGHT 

The critical feature of this application is the analysis of 
the inviscid flow, since known two-dimensional boundary- 
layer solutions can be readily employed once this flow is 
known. Accordingly, the following discussion is restricted 
to the inyiscid flow problem. 

Now it was shown previously in this paper that a large 
class of hypersonic flows which are basically three-dimen- 
sional can be calculated with a generalized shock-expansion 
method which is analogous to that employed in reference 7 
for studying flow about airfoils. Speciflcally, this treatment 
is permissible when disturbances associated with the diver- 
gence of streamlines in planes tangent to a surface can be 
considered negligible compared to those associated with the 
curvature of streamlines in planes normal to the surface (see 
eq. (26)). For the case of nonindined bodies of revolution 
which are curved in the stream direction, this requirement is 
satisfied when the hypersonic similarity parameter K is 
greater than about 1 (see ref. 10). For inclined bodies, an 
additional restriction is imposed. This point is perhaps best 
clarified by considering the problem of calculating flow at 
the surface. 

FLOW AT THE SURFACE 

It follows from the inviscid flow analysis that when the 
generalized shock-expansion method applies in the region 
downstream of the vertex, surface streamlines can be 
approximated by geodesic lines. The only geodesics on 
the surface of a body of revolution which, like streamlines, 

413072 — C7 70 


do not intersect each other are the meridian lines. In 
addition, the meridian lines are the only geodesics which, 
like the streamlines, pass through the vertex. When the 
shock-expansion method is applied, then, surface streamlmes 
are approximated by meridian lines. Strictly speaking, 
however, this approximation is valid only in the case of 
a:«l. (It is, of course, always true, independent of a 
on the extreme windward and leeward sides of a body.) 
Evidently, then, the generalized shock-expansion method 
should be applicable to curved bodies of revolution only at 
small angles of attack in flows characterized by a value of 
the hypersonic similarity parameter greater than about 1. 

The procedure for determining flow conditions at the 
surface of a lifting body is entirely analogous to that em- 
ployed in the application of the shock-expansion method 
to the nonliftmg body (ref. 10). Thus, it is assumed that 
the flow at the vertex is the same as that for a cone tangent 
to the body at this point and, hence, may be determined 
from existing conical-flow theory (see, e. g., ref. 31 for 
moderate supersonic Mach numbers and ref. 32 for high 
supersonic Mach numbers). More specifically, the Mach 
number at the vertex under the vortical layer “ may be 
calculated by means of the pertinent conical-flow expressions 
in reference 31 or reference 32. The variation of Mach 
number downstream of the vertex is then obtained by means 
of the Prandtl-Meyer angle v (see, e. g., ref. 33) which in 
turn is determined from the isentropic expansion relation 

5a+»’a=5b-J->'jj (41) 

where A and B are different points on the same meridian 
line. Since the flow is isentropic in the windward plane of 
symmetry downstream of the shock at the vertex and around 
the surface of the body, the pressure distribution (in co- 
eflficient form) is readily obtainable with the aid of the expres- 
sion 



where p./p,, is the pressure rise across the shock at the 
vertex on the windward side of the body and is determined 
from conical-flow theory. The ratio p/p, is given by 


JL. 

P. 




y 

7-1 




(43) 


where M is known from equation (41) and M, is the Mach 
nmnber unmediately downstream of the shock on the wind- 
ward side of the body at the vertex and, hence, is also 
determined from conical-flow theory. 


FLOW OFF THE SURFACE 

Flow in meridian planes around bodies of revolution may 
be calculated by the generalized shock-expansion method in 
much the same manner as the procedure employed in refer- 
ence 7 for flow about airfoils. However, the application of 
the method is somewhat more complicated for the case of a 


» since a vortical layer exists around the body gnrlace at the vertei (see, e. g., ret 13) a 
vortical layer must also exist downstream oT the vertex. 
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body of revolution since now the influence of the conical- 
flow region at the vertex must he considered. An analysis 
for flow in the region of the vertex of a nonlifting body of 
revolution (AC>1) presented in reference 34 and expres- 
sions were developed which yield the shock-wave curvature 
as well as flow conditions along a line (normal to the hody 
axis) a short distance downstream of the vertex. This 
analysis was extended to lifting bodies and more general 
expressions were presented in reference 32. Thus, initial 
conditions m the region of the vertex can he established. 
There remains the determination of flow conditions along 
meridian Imes downstream of the vertex externally adjacent 
to the vortical layer. These conditions may be determined 
in the same manner as for flow directly on the surface (i. e., 
under the vortical layer), except that now initial flow condi- 
tions externally adjacent to the vortical layer at the vertex 
are employed in the isentropic expansion relations. Con- 
struction of the flow field between the shock and the vortical 
layer in each meridian plane can then proceed in a manner 
analogous to that for the two-dimensional airfoil discussed 
in reference 7. To fllustrate, consider the flow in a meridian 
plane of a lifting body of revolution (see fig. 6). 

y 



Figubb 6. — Schematic diagram of flow field about a body of revolution. 


FORCES AND MOMENTS ACTING ON BODIES OF KEVOLUTION 

It is of interest now to consider briefly the forces acting 
on a hody of revolution. In the previous discussion, 
attention was called to the fact that the flow is isentropic 
in the windward plane of symmetry at the vertex as woU as 
on the surface of the hody. This result materially reduces 
the net labor associated with carrying out the calculations 
to determine the pressure distributions around the body 


downstream of the vertex since the pressure rise, as well as 
the change in entropy through the shock, need he considered 
only in this plane at the vertex. The normal-force, axial- 

force, and pitching-moment coefficients may he 
from the expressions 

obtained 


(44) 

X * (X 0 * 

(46) 

16 r* V 

(46) 


respectively, where d is the diameter of the base, r is the 
radius of the hody, ip is the meridian angle measured from 
the plane of symmetry on the windward side of the body, 
and X is measured along the body axis. If equations (44) 
and (46) are differentiated with respect to a and the condi- 
tion of constant entropy on the surface is employed, there 
is obtained 


fp \ fT' ^ 

V Jo Vtr sin2/n &a 

and 



r cos (p dtp dx 
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AH fluid properties at points N, A', A, B and so forth, on 
the body surface external to the vortical layer are calculated 
with the aid of the obUque shock-wave, conical-flow, and 
expansion equations. Mow conditions along the line AC 
may then he determined (see ref. 32). It will be recalled 
(see ref. 7) that a basic condition employed in constructing 
flow fields about airfoils by the generalized shock-expansion 
method is that the pressure is constant along Mach lines 
emanating from the surface. In the case of flow about 
pointed bodies of revolution, this condition can be relaxed 
to account for the small variations in pressure due to the 
influence of the conical type flow in the region of the vertex. 
The procedure is as follows. The Mach line A'C is con- 
structed using the known conditions in the region NAG 
shown in the sketch. The net pressure change along this 
Mach line (i. e., Po~ 3 >a') is thus determined. This pressure 
difference is then assumed to represent the net pressure 
change between the hody surface and the shock along each 
Mach line emanating from the surface downstream of the 
vertex. The flow field is constructed using this criteria in 
conjimction with the isentropic expansion relations for flow 
along stream lines. 


'n \ 16 f' p sm 2 /i^, d /p/f\ , 7 

J = Ar a- 3 i !7 — • o ST" ( ~ ) rx cos <p dtp dx 
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These expressions define the initial slopes of the normal- 
force- and pitching-moment-coefficient curves, respectively, 
and may be rewritten in terms of the initial nonnal-force- 
curve slope for a cone tangent to the body at the vertex; 
thus. 


(i) ^ 


and 




sin 2 fty 
sin 2 /i 






tan 



p sin 2 /X 7 V 
P}f sin 



where the subscript TON refers to a cone tangent at the 
vertex. The calculations necessary to determine the initial 
normal-force- and pitching-moment-curve slopes for a body 
of revolution are then relatively simple, since ipifc)ro!f 
be easily obtained from reference 12 or from chart 8 in 
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reference 33. The Mach numher and pressure distribu- 
tions along the body are obtained by the shock-expansion 
method for the case a=0. When these distributions have 
been determined, the integral terms in equations (49) and 
(60) are easUy evaluated by numerical integration or by 
graphical methods. 

SIMPLIFIED BXPHESSIONS FOR SLENDER BODIES 

In the case of slender bodies traveling at very high super- 
sonic Mach numbers and very small angles of attack the 
calculations of fluid properties at the surface become rela- 
tively simple. In fact, fluid properties downstream of the 
vertex may be related to those at the vertex by means of 
explicit algebraic expressions. In particular, the local Mach 
number and pressure distributions on the surface of a slender 
body may be written (see ref. 32) 





respectively, where 5 is measured relative to the body axis 
in these and subsequent expressions. Equations (61) and 
(62) combine with the corresponding conical-flow equations 
(ref. 32) to predict the ratios of local to free-stream Mach 
numbers and local to free-stream static pressures to be the 
same at corresponding points on related bodies, provided the 
flow flelds about these bodies are defined by the same re- 
spective values of the hypersonic similarity parameters 
and M^a (or a/Sn). These predictions are in agreement 


with those of reference 3 for inviscid flow about slender 
three-dimensional shapes, and they enable the solution of 
equations (61) and (62) in terms of tabulated functions of 
the simil arity parameters. Calculations over a range of 
McS^r from 0.60 to “ and a/8x from 0 to 1 were carried out 
for flow at the vertex of a body of revolution*® and the results 


of these calculations 
(P./PJ^. 


for the flow parameters 


(MX 


Kv 


(MM 


0 


and (M,5jy)^_o are tabulated m table I for 30° 


increments of p from 0 to ir. For a given Ma,di^ and 
Ma.ct, the Mach number on the surface of a body downstream 
of the vertex is readily obtained with the aid of these tabu- 
lated parameters when used in conjunction with equation 
(51). The pressure coefficient is easily calculated by means 
of equation (52). 

The results from table I may also be used to good advan- 
tage in determining the initial slopes of the normal-force- and 
pitching-moment-coeffident curves for slender bodies.*^ For 
example, when AC>>! and 5<<1, equations (49) and (50) 
combine with equations (51) and (52) to yield 


(Otr^a-o — 1^1 ^ (AfjiAr) ^1 — 
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respectively, where 
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1+y^ (Mj^r 
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(55) 


These expressions are easily evaluated with the aid of the 
tabulated flow parameters in table I for the case a/Sif=0. 

EXPEREMENT 

In order to obtain a check on the predictions of the pre- 
ceding theoretical analysis, the pressures acting on the 
surfaces of bodies of revolution corresponding to values of 
the hypersonic similarity parameter K from 0.60 to 2.1 at 
Mach numbers from 3.00 to 6.30 were determined experi- 
mentally, The bodies were tested at angles of attack up 
to 16°. A brief description of these tests follows. 

TEST APPARATUS 

Tests were conducted in the Ames 10- by 14-inch super- 
sonic wind tunnel. A detailed description of the wind tunnel 
and auxiliary equipment may be found in reference 35. 


The pressures acting on the model surfaces were measured 
with a mercury U-tube manometer or by means of McLeod 
gages when the pressures were low enough to be recorded 
on the latter. 

Pressure-distribution models were mounted on a 0° model 
support and on 5°, 10°, and 15° bent supports. The test 
models were two tangent ogives having fineness ratios 3 and 
5 and two cones having the same vertex angles as the ogives. 
The dimensions of these models and location of the pressure 
orifices are shown in figure 7. 

The conical flow expressions presented in reference 32 were employed in these calculations. 

It win be noted in table I that the value of not given for all values of p. This 

results from the fact that as ar-*iN and jVf ® i/r-* ® , the assumption of an infinitesimally thin 
vortical layer is violated and, hence, the slender-oone theory yields onrealistlo results for these 
oondJtlons. 

w The initial axlal-fbrce-cnrvo slope Is, of course, icro due to symmetry. 
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Orifice location shown by + 
All dimensions in Inches 



Circular- arc ogival models Conicol models 

Figttbh 7. — Dimensions of pressure distribution test models showing location of pressure orifices. 


TESTS AND PEOCEDDKE 

Pressures on the model surfaces were meastired at 0°, 
6°, 10°, and 15° angles of attack and at test Mach numbers 
of 3.00, 4.25, and 5.05. Pressures on the fineness-ratio-3 
ogive (as well as on the corre^onding cone) were also obtained 
at a test Mach niunber of 6.30 and at 0° and 5° angles of 
attack. The Reynolds mnnbers (based on maximum diam- 
eter of the ogives) were 1.09 million at Mach inmabers 3.00 
and 4.25, 0.52 million at Mach nmnber 5.05, and 0.22 
million at Mach number 6.30. 

The pressures arovmd the cone surface (0° to 360°) at 
meridian stations 45° apart were recorded simultaneously 
at each Mach number and angle of attack. In the case of 
the two ogival models, the pressures were recorded at 
meridian stations 90° apart. Each model was then rotated 
45° about its longitudinal axis (except at 0° angle of attack) 
and the process repeated. 

ACCTJRACrr OF TEST RESULTS 

In the r^on of the test section where the models were 
located, the variation in Mach nmnber did not exceed ±0.02 
at Mach numbers from 3.00 to 5.05 and ±0.04 at Mach’ 
number 6.30. 

The precision of the computed pressure coefldcients was 
affected by inaccuracies in the pressure measurements, as 
well as uncertainties in the stream angle and the free^stream 


dynamic pressure. The restilting errors in the pressure 
coefficients were generally less than ±0.005 throughout 
the Mach number range for aU angles of attack.^ 

COMPARISON OF THEORY WITH EXPERIMENT AND DIS- 
CUSSION OF RESULTS 

According to inviscid theory, the hypersonic similarity 
parameter, 5T, is a significant index to when the generalized 
shock-expansion method can be used to calculate three- 
dimensional flows. It was indicated in this connection that 
the generalized method should be applicable to bodies of 
revolution when K is greater than about 1. In order to 
check this prediction, the pressure distributions on the 
surfaces of two ogives (having fineness ratios 3 and 5) 
traveling at Mach numbers 3.00, and 4.25, and 5.05 and 
at angles of attack of 0°, 5°, 10°, and 15° were calculated 
by the methods of this paper. Pressure distributions on 
the fineness-ratio-3 ogive at angles of attack of 0° and 5° 
were also calculated for a Mach nmnber of 6.30. The 
conical-flow theory presented in reference 32 was employed 
in these calculations for determining initial flow conditions 
at the vertices of both the lifting and nonlifting bodies. 

Comparing first the predictions of theory with experiment 

u T hB experimental data shown in figure 8 and aU sabsequent figures represent on ovcrogo 
of the pressures recorded at each station on a body. The scatter of data was Inconsequentially 
nrtflTi (seereL 32). 
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(a) AT„=3.00 (b) iVI.=4.26 (c) Af„=5.06 (d) M„=6.30 

FiGTmB 8. — Variation of pressure coefficient along ogives at a=0°. 


for the case of zero lift, we obserye in figure 8 that the shock- 
o-xpansion method predicts surface pressure coefficients close 
to those obtained experimentally at values of K greater than 
1. As would be expected, too, the agreement between the 
predicted coefficients and experiment tends to improve up to a 
Mach number of 5.05. The results of a characteristics solu- 
tion for a fineness-ratio-3 ogive at Af„=3.00 (from ref. 36) 
are also shown for comparative pm^poses. Characteristics 
solutions are not available for the other cases; however, the 
results of Rossow (obtained by correlating the pressures 
yielded by characteristics solutions according to the hyper- 
sonic sirml.arity law; see ref. 37) are shown. It is evident in 
this figure that the agreement between these results and 
those 3 uelded by the shock-expansion method improves with 
increasing K over the Mach number range presented. At the 
highest Mach number of 6.30 we observe, however; that both 
methods yield pressure coefficients which, although in agree- 
ment, are appreciably lower than experiment. There is no 
particular reason, on the basis of past experience or otherwise, 
to doubt the accuracy of the characteristics theory for this 
body. In this connection, it should be noted that the theory 
is generally in good agreement with experiment at all the 
lower Mach numbers. It seems logical, therefore, to suspect 
that the departure of theory from experiment at Ma>=6.30 is 
caused by viscous effects in the flow. More specifically, it ie 
suggested that this departure may be traced to a substantial 
increase in thickness of the laminar boimdary layer on the 
ogive. The low Reynolds number of the tests and, to .a 
somewhat lesser extent, the high Mach nmnber could produce 


such an increase. This matter will be considered further in 
the discussion of hypersonic boundary-layer calculations pre- 
sented later in the paper. 

It is appropriate now to consider the reliability of the 
shock-expansion method for lifting bodies. As shown in 
figures 9, 10, and 11, the theory yields good agreement with 
experiment on the windward side of the fineness-ratio-6 ogive 
except at Af„=3.00 (£^=0.60).'® Disagreement is evident, 
however, on the leeward side of the body at all Mach num- 
bers. In the case of the fineness-ratio-3 ogive (figs. 12, 13, 
and 14), agreement is generally better over the entire body 
at each angle of attack, particularly at the higher values of K. 
It will be recalled from figure 8 that at a=0° the longitudinal 
pressme distributions on both ogives indicated that the accu- 
racy of the shock-expansion method increased as increased. 
Figures 9 through 14 indicate that, as would be expected, this 
trend carries over to the case of lifting bodies. It is inter- 
esting to note, also, that reasonably good agreement with ex- 
periment is obtained when FC>1 even though Ac- 

cordingly, it is suggested that so long as .K!> 1 and a/S^<C^l, 
the generalized shock-expansion method can be employed .to 
predict surface pressures along meridian lines as though they 
were streamlines, with little sacrifice in accuracy. In this 
connection, it should be noted that the meridian lines on the 
extreme windward and leeward sides of the body (i. e., <p=0° 
and ip= 180°, respectively) are exactly streamlines. 

“ tt xhoiild be noted in figure 9 that Stone’s second-erder sohitlon is employed at the vertoi 
since the oonlcal-fiow theory of reference 3J is not applicable for these conditions (J, d, 
Af CO -3.00 and Jw- 11.43°). 
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It appears in figures 9 through. 14, that the most important 
factor influencing the accuracy of the shock-expansion 
method is the reliability of the conical-floiv theory, since the 
inaccuracies at the vertex appear to be reflected strongly in 
the pressures downstream of the vertex. The question 
naturally arises, then, how good are the predictions of the 
method when experimentally determined initial conditions 
at the vertex are employed? To answer this question, the 
pressure coefficients on the surfaces of the two ogives imder 
discussion were determined in the following manner. Initial 
conditions at the vertex were determined from the measmed 
static pressures around a cone (corresponding to the vertex 
angle of the body) in conjimction with the measured shock- 
wave angle (in the plane ^s>=0°) obtained from schheren 
photographs of the conical flow field. The pressure co- 
efficients downstream of the vertex were then calculated as 
before. The results of these calculations for Mach numbers 
3.00, 4.26, and 6.05 are compared with experiment in figirres 
16 and 16 for a=16°. Results for a=15° are presented 
because at this angle of attack the applicability of the 
conical-flow solutions is most mai^inal. It is observed in 


figure 15 (a) that the theory yields results which indicate an 
underexpansion of the flow on the sides of the body (vj=46° 
and 95=90°). This result is not surprising since a/5«-)>l 
and K<^1. It would be expected, then, that the true stream- 
lines would deviate considerably from a meridian line. In 
other words, flow disturbances in planes tangent to the body 
at the surface are no longer small compared to those in axial 
planes. It can be seen from figures 15 (b) and (c) that as 
the Mach number, and hence, K, is increased, bettor agree- 
ment is obtained. This result is attributed in part to the 
fact that the streamlines of the flow deviate less from meridian 
lines as is increased. The same general trend may be 
noted in figure 16 for the fineness-ratio-3 ogive. However, 
in this case, and over-all agreement between theory 

and experiment is improved. In fact, good results are con- 
sistently obtained by theory except on the extreme leeward 
side of the body where it is probable that viscous effects are 
influencing the pressures. There may be some separation 
of flow over this portion of the body although no evidence of 
this could be determined from the schlieren photographs. 
In the case of the fineness-ratio-5 ogive, schlieren evidence 
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indicated flow separation on the leeward side of the body for 
all Mach numbers at «=16°. It is evident from these 
figures that in any event the shock-expansion method will 
yield better results when initial conditions at the vertex are 
determined from cone tests rather than from presently 
available cone theory. 

There now remains the determination of the accuracy of 
the predictions of the generalized shock-expansion method 
for the flow field (other than the siuface) about a lifting 
body of revolution. To this end, fiow in the plane of sym- 
metry (^=0° and ¥>=180°) was calculated for each ogive 
traveling at a Mach number 6.05 and at an angle of attack 
of 10°. Flow in a side meridian plane (v5=90°) was also 
calculated for the fineness-ratio-3 ogive. The resulting 
shock-wave shapes are compared with the actual shapes (ob- 
tained from scblieren photographs) in figure 17. The theo- 
retically determined conical shocks are also shown for con- 
trast. In the case of the fineness-ratio-3 ogive (iir=1.68 and 
a/Sff=0.5Z), theory and experiment are observed to be in 



Meridian 


excellent agreement in the plane of sy mm etry, The same 
observations may be made for the side meridian plane. In 
this latter connection, it is of interest to point out that es- 
sentially the same result is obtained when the shock is as- 
sumed circular in cross-sectional planes and its location de- 
termined from the calculations in the plane of symmetry. 
In the case of the fineness-ratio-5 ogive, the poor agreement 
on the leeward side of the body is due to the limitations of 
the conical fiow theory employed at the vertex. If experi- 
mentally determined initial conditions are employed good 
agreement with experiment downstream of the vertex is 
obtained. 

Although the predictions of the generalized shock-expan- 
sion method have been checked only at the inner and outer 
boundaries of the fiow field, it is expected that equally good 
results would be obtained at intermediate points in the fiow 
field. This conclusion is based on the fact that the bow 
shock waves were obtained as a result of the calciilations of 
these intermediate points. 



angle, <j>, deg 


Fioubb 10. — Circumferential variation of pressure coefficient on a fineness ratio 6 ogive at Af„=4.26; K=0.85, 





Pressure coefficient, Cp 





Pressure coefficient, Cf. 
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FiauBB 12. — Circumferential variation of pressure coefficient on a fineness ratio 3 ogive at ikf„=3.00; fiT=1.00. 



Pressure coefficient , 
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Meridcon angle, <j>, deg 


iPiQUEB 13. — Cirovunferential variation of pressure ooeflSoient on a fineness ratio 3 ogive at Af„=4.26; If ■=1.42, 



Pressure coefficieef, Cp 
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30 60 90 120 ISO 180 0 30 60 90 120 ISO 180 


Meridian angle, deg 

Piatma 14, — Circumferential variation of pressure coefficient on a fineness ratio 3 ogive at Af„=6.06; ,£^=1.68. 
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Symbol x/t 

o 

0 (Cone doto) 

□ 

.27 


54 

’A 

50 


1.00 


Shock-expansion method 
employing experimentally 
determined cone conditions 


30 60 90 120 150 180 0 30 60 90 120 150 180 O 30 60 90 120 

Meridian angle, (f>, deg 

(a)if„=3.00;S:=0.60 (b) ilf„=4.26;X=0.85 (c) =6.06; S= 1.01 

Fioxjbb 16. — Circumferential variation of pressure coefl5cient on a fineness ratio 6 ogive at a=16°. 


Symbol x/1 
o 0 

a .22 

O .44 

a .66 

a 1.00 


Shock-expansion method 

^ employing experimentally 

\ I determined cone conditions 


Meridion angle, deg 

(a) ilf„=3.00; Kr=l.QQ (b) Af„=4.26; S:=1.42 (c) Af„=6.06; Z=1.68 

FiatTEB 16. — Circumferential variation of pressure coefficient on a fineness ratio 3 ogive at a=16°. 
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Experimental nose solution 





Generalized shock-expansion 

Theory method 

Conical shock 

Experiment o Determined- from schlieren 
photographs 


(b) K=1.01 


Fioube 17. — Shook-wave shapes for two ogives at Af„=6.06 and a=10°. 
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It is appropriate now to consider briefly the forces exper- 
ienced by the ogives. To this end, normal-force coefficients 
were obtiiined by integrating the theoretical pressure distri- 
butions for the two ogives at a Mach number of 6.05. The 
results of these calcrflations are compared in figure 18 with 
those obtained from integrated experimental pressure dis- 
tributions for values of K. of 1.01 and 1.68. It is observed 
that although theory yields results which are, in general, 
higher than those obtained by experiment, agreement im- 
proves with increasing K. The same trend with K is evi- 
dent for the initial normal-force-curve slopes obtained with 
the aid of equation (49). Axial forces were also obtained 
for these ogives and, as indicated in figure 18, the shock- 
expansion method yields generally good agreem en t with 
experiment even at a value of as low as 1. 



(a) Z=1.01 

(b) S:=1.68 

FiGtTBB 18. — ^Normal- and axial-force coeflBcIents for ogives at 

M„=6.05. 


efficients for the two ogives at a=0° and a=6° are com- 
pared with experiment in figures 19 and 20. It appears from 
a comparison of figures 8 and 19 that the slender-body theory 
will yield more accurate drag coefficients than the general 
theory at a=0°, particularly at the lower values of K. 
This result is, of course, fortuitous. In the case of lifting 
bodies (fig. 20) the slender-body theory yields results which 
are somewhat less satisfactory at aU values of K, However, 
the theory displays sufficient accuracy for many engineering 
purposes even at K=l. ' This point is particularly evident 
for the more slender of the two bodies as indicated in figure 
20. It is also interesting to note the comparison of theory 
and experiment shown in figure 21 for the initial normal- 
force-cmwe slopes and centers of pressme of a family of 
ogives at Mach numbers from 3.00 to 6.30. The experi- 
mental data were obtained in the Ames 10- by 14-inch 
supersonic ■wind tunnel. There is good correlation of these 
data ■with Mod^, the hypersonic similarity parameter for 
slender bodies, and there is good agreement ■with theory for 
values of greater than 1. In view of its simplicity, 

then, the hypersonic slender-body theory should prove useful 
and its application is further facilitated by the presentation 
in this paper of tabulated values of the pertinent flow para- 
meters for selected values of M^Sjf and a/Sx (see table I). 

Up to this point we have been concerned almost entirely 
■with the in^viscid theory and its comparison "with experi- 
mental data relatively free of effects of viscosity. As a final 
point, it is appropriate to test the two-dimensional boundary- 
layer concept of this paper by considering flows which are 
s ignifi cantly influenced by viscous effects. 

In this connection it was noted early in the pre^vdous dis- 
cussion that in'viscid theory yielded pressure coefficients 
which were substantially lower than experiment at M„= 
6.30. This discrepancy was traced to the thick laminar 
boundary layer on the test body. According to theory it 
should be possible to calculate this boundary layer approxi- 
mately by means of simple two-dimensional techniques. 
This possibility was checked by calculating the laminar 
boundary layer on the fineness-ratio-3 ogive at a Mach 
number of 6.30 and angles of attack of 0° and 6°. The 
two-dimensional theory of reference 38 was employed.’^ 
The body ordinates were increased by an ammmt equal to 
the displacement thickness of the boundary layer. The 
pressure distribution about the distorted body was then 
obtained ■with the generalized shock-expansion method. 
These corrected pressure distributions and the original 
uncorrected distributions are presented in figure 22 along 
■with experiment. It is observed that Avhile imcorrected 
press'ures are definitely low, the corrected pressure distribu- 
tions are in good agreement ■with experiment. In this case, 
then, relatively simple two-dimensional methods of cor- 
recting pressure distributions for the presence of the 
bo'undary layer are, as indicated by theory, applicable to 
the body of revolution. 


Let us consider now the predictions of the hypersonic 
slender-body theory. To this end, calculated pressure co- 


« The theory breaks down at the vertex of the body, much as at the leading edge of an air- 
fofl. It Is therefore not applied In thl< region, and, consistent with a praotlco successfully 
employed with airfoils, viscous effects are Ignored In cal cola ting flow at the vertex. 





Pressure coefficient, 
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Figubb 19. — Pressure distributions predicted by hypersonic slender-body theory for ogives at a=0°. 





Pressure coefficient, 




Meridian ongle, <t>, deg 


(a) M„=3.00 M„=6.05 

Figubb 20— Pressure distributions predicted by hypersonic slender-body theory for ogives at «=6' 



FiGimB 21. — Normal-force derivatives and centers of pressure for ogives 
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(a) a=0° 

(b) a=5° 

PioiniE 22. — Effect of boundary layer on pressures acting on a fineness 
ratio 3 ogive at Af„=6.30. 

CONCLUDING REMARKS 

A method of characteristics employmg pressure and flow 
inclination angles as dependent variables was used to 
obtain a simplified approximate method for calculating 
three-dimensional flows at high supersonic speeds. It was 


fmmd that when the flight Mach nmnber is sufficiently 
large compared to 1, flow in the osculating planes of stream- 
lines in regions free of shock waves may be of the generalized 
Prandtl-Meyer type — surface streamlines in this event may 
be treated as geodesics. In the case of slender bodies, these 
results apply to nonsteady as well as steady flows. The 
two-dimensional approach to three-dimensional hypersonic 
flows was also extended to steady boundaiy-layer flows. 

Bodies of revolution in steady hypersonic flight were 
considered as an example of shapes producing three-dimen- 
sional flow fields which appear locally two-dimensional. 
With the assumption of conical flow at the vertex and small 
angles of attack, simple approximate solutions were obtained 
which yield the Mach number and pressure distributions on 
the surfaces of such bodies. Surface streamlines were 
approximated by meridian lines and the flow field in meridian 
planes was calculated by means of a generalized shock- 
expansion method. In the special case of slender bodies, 
simple explicit expressions were obtained for the Mach 
number and pressure distribution on the surface. 

Surface pressures and shock-wave shapes were obtained 
experimentally at Mach numbers from 3.00 to 5.05 for two 
ogives having fineness ratios 3 and ■ 5 and for two cones 
having the same vertex angles as the ogives. The predic- 
tions of the methods of this paper for the surface pressures 
and shock-wave shapes were found to be in good agreement 
with experiment at values of K of about 1, or greater, when 
albif (the ratio of angle of attack to semivertex angle) was 
about 1/2 or less. For increasing values of this parameter, 
agreement deteriorated but was still reasonably good for 
values of ajS^ up to about 1. Experimental surface pres- 
sures at a Mach number of 6.30 and angles of attack of 0° 
and 5° were also obtained for the fineness-ratio-3 ogive. 
The predictions of the shock-expansion method when 
employed in conjunction with a two-dimensional boundary- 
layer calciflation were found to be in good agreement with 
experiment. 

In view of these results, it is concluded that the generalized 
shock-expansion method should prove usefffi in treating 
three-dimensional hypersonic flow fields about practical 
aerodynamic configurations. Fmthermore, it is indicated 
that methods of treating two-dimensional hypersonic 
boundary layers may, in like manner, prove useful in pre- 
dicting three-dimensional hypersonic boimdaiy layers. 

Ames Aeronautical Laboratory 

National AnvrsoRT Committbe for Aeronautics , ■ 
Moffett Field, Calif., Aug. 16, 1952 
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TABLE L— TABLE OP FUNCTIONS FOR HYPERSONIC SLENDER-BODY METHOD 











































































